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Abstract 

We investigate the breakdown of supersymmetry at finite temperature. 
While it has been proven that temperature always breaks supersymmetry, 
the nature of this breaking is less clear. On the one hand, a study of the 
Ward-Takahashi identities suggests a spontaneous breakdown of supersym- 
metry without the existence of a Goldstino, while on the other hand it has 
been shown that in any super symmetric plasma there should exist a mass- 
less fermionic collective excitation, the phonino. Aim of this work is to 
unify these two approaches. For the Wess-Zumino model, it is shown that 
the phonino exists and contributes to the supersymmetric Ward-Takahashi 
identities in the right way displaying that supersymmetry is broken sponta- 
neously with the phonino as the Goldstone fermion. 
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1 Introduction 



Supersymmetry has become a central element in the extensions of the standard 
model of elementary particle physics that currently attract most attention. It 
is not only an intriguing mathematical concept with a distinct importance for 
theoretical physics as the only nontrivial extension of the Poincare algebra in 
relativistic quantum field theory, but it has some very attractive properties also 
from the point of view of phenomenology 

On the other hand, our insight into the evolution of the universe has made 
an enormous progress over the last decades, and there has been a stimulating 
interchange of ideas between cosmology and quantum field theory. Of particular 
importance is the behaviour of symmetries in a hot plasma like it is believed to 
have existed shortly after the big bang, since phase transitions may have left their 
traces in the present state of the universe. For example, it is well-known that 
the spontaneously broken gauge symmetry of the standard model is restored at 
high temperature, a phenomenon that appears quite generally for global and gauge 
symmetries. So, if we believe that supersymmetry was realized in the hot early 
stages of the universe in one or the other way, it is indispensable to understand its 
behaviour at high temperatures. 

Considering the fact that this is such a fundamental problem, it is astonishing 
how little has been done to solve it. Although it was realized that supersymmetry 
generally breaks down at finite temperature already in the earliest works on the 
subject PQI^l, the subsequent literature contains many contradictory statements, 
and our understanding of the mechanism of supersymmetry breaking at finite 
temperature is still unsatisfactory. 

The current status can be summarized as follows. Only recently, it has been rig- 
orously proven that supersymmetry is always broken at any finite temperature j^j. 
In fact, this is not surprising, as in thermal field theory the ground state is described 
by a statistical ensemble with different populations of bosons and fermions. Since 
the thermal ground state is responsible for the breakdown of supersymmetry, it has 
much in common with a spontaneous breaking. Thus, a natural question to ask is 
whether it is associated with the existence of a massless Goldstone fermion. The 
investigation of the supersymmetric Ward-Takahashi identities j3J [3] showed that 
there must be a zero-energy Goldstone mode. However, since the rest frame of the 
heat bath also breaks Lorentz invariance, this mode is not necessarily associated 
with a propagating Goldstone particle. Nevertheless, in jB] (clarifying earlier ideas 
in [7]) it has been shown from a complementary point of view that any model with 
thermally broken supersymmetry should contain a massless fermionic collective 
excitation, similar to the appearance of sound waves in a medium with sponta- 
neously broken Lorentz invariance. The associated particle was baptized phonino 
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because of its similarity to the phonon. For the simple case of the Wess-Zumino 
model, the existence of this phonino was proved. 

In this work, we will again focus on the Wess-Zumino model. Though it ap- 
pears relatively simple and has already been investigated in a number of papers, 
it reveals a very interesting structure. Our aim is to unify the different results in 
the literature and to obtain a complete picture of the breakdown of supersymme- 
try in this model. To this end, we will give an explicit proof of the existence of 
the phonino and investigate its contributions to the Ward-Takahashi identities of 
broken supersymmetry. It turns out that supersymmetry is indeed broken spon- 
taneously by the heat bath with the phonino playing the role of the Goldstone 
particle. The results obtained for the Wess-Zumino model allow us to infer the 
behaviour of more general models. 

The paper is organized as follows. After briefly reviewing the theoretical frame- 
work and the Wess-Zumino model in sections El and HU the established knowledge 
about supersymmetry and one-loop behaviour of the Wess-Zumino model will be 
presented and partly extended for our purposes in sections E| and 03 In section El 
we present a full calculation of the fermion propagator for low momenta and give 
an explicit proof of the existence of the phonino. Finally, section [7| is devoted to 
the investigation of the Ward-Takahashi identities. 



2 Thermal field theory 

We will work in the framework of finite temperature field theory which is the 
appropriate formalism for the description of quantum fields in thermal equilibrium. 
The thermal background will be described by the canonical ensemble with a density 
matrix 

p = Z-'e-^. 

Here, H is the Hamiltonian, and the inverse temperature (3 = (k^T)' 1 is chosen in 
units so that Boltzmann's constant is unity. The partition function Z normalizes 
the density operator so that the expectation value of an observable O is given by 

{0) p = \xpO. 

The basic effect of the presence of the thermal background on free quantum 
fields is a modification of the propagators. Since annihilation operators do not 
annihilate the thermal ground state, the propagator of a scalar field A gets an 
additional thermal contribution. It reads 

D(p)= fd A xe^ x ~ y \TA{x)A{y))p = - + 2vr 6(p 2 - m 2 ) n B (po), (1) 
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where 



™b(Po) 
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(2) 



e^bol _ i 

is the Bose-Einstein distribution function. Analogously, the free fermion propaga- 
tor is given by 



S(p)= U A xe ip{x ~ y) (T^(x)^(y)) p 



i{j/> + to) 

p2 _ m 2 _|_ l e 

with the Fermi-Dirac distribution function 

n F (Po) 



1 



e^bol + 1 ' 



' + to) 2n5{p 2 — m 2 ) fi-pipo), 

(3) 

(4) 



Already at this point it becomes clear that supersymmetry has a hard time at 
finite temperature since the different distribution functions lead to quite different 
thermal contributions to the bosonic and fermionic propagators. 

The treatment of an interacting theory requires much more effort. From the two 
formalisms that have been developed for this purpose, we will make use of the real- 
time formalism since it is suited for the direct perturbative calculation of thermal 
correlation functions in real time by standard Feynman diagram techniques. For 
consistency, this formalism requires a doubling of the degrees of freedom. Formally, 
one introduces for each field a ghost field with the same interaction as the original 
field, only of opposite sign. The thermal propagator of a scalar field is then given 
by a matrix of the form 



D{p) = M(p 



where 



M(po) 



/ cosh#(po) 
I sinh 9(po) 



p 2 — m 2 + ie 




sinh 9(po) 
cosh 9(po] 





i 



\ 



M(po), 



(5) 



' 2 — m 2 — ie / 



V 



with sinh^fjoo) = ^b(Po 



The 11-component of this matrix gives back the thermal propagator ((TJ. More 
general correlation functions can be calculated perturbatively by evaluating the 
same diagrams as in the vacuum theory, where external lines always correspond to 
physical 1-fields. It will later turn out that the ghost fields do not lead to relevant 
contributions in our calculations, but neverless must be considered in a consistent 
treatment. 

The matrix structure (j3J) allows also the understanding of self energy correc- 
tions. Without going into any details (for which we refer to the standard literature 
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as jH]), we note that the full propagator must have a similar structure, and the 
thermal self energy corrections only lead to the usual shift 

p 2 — m 2 — > p 2 — m 2 — Hp(p), 

while the overall structure (JIJ of the thermal propagator remains unchanged. The 
real part of the thermal self energy function Hp(p) is directly accessible to a per- 
turbative calculation. It coincides with the real part of the self energy function 
with physical external lines, while the imaginary part requires minor corrections. 

In case the self energy is real and small compared to the mass, it only induces 
a small, temperature-dependent mass shift. One should however note that there 
is now an ambiguity in defining the mass. Because of the breakdown of Lorentz 
invariance in the heat bath, the self energy Hp(p) not only depends on p 2 but also 
on p. Therefore, the value of the propagator at zero momentum or its pole for 
vanishing three-momentum lead to different notions of mass while they coincide 
at zero temperature. For reasons that will become clear later, we will adopt the 
definition of mass as the pole of the full propagator in the limit of vanishing three- 
momentum, 

m 2 p = m 2 + Hp(m, 0). (6) 

The fermionic case appears quite similar. Here, the thermal self energy leads 
to a temperature-dependent shift 

$ -m^i> -m- Ep(p) 

so that the thermal mass is given by 

m 2 = ^tr[(l + 1 °)Zp(m,0)], (7) 

as long as the self energy is small. In general, the propagator can of course have a 
much more complicated structure and in particular involve an imaginary damping 
part. 

In finite temperature field theory, ultraviolet divergences can be absorbed by 
the same redefinition of the parameters as in the vacuum theory so that higher 
order vacuum contributions are small after renormalization. In this work, we will 
only be concerned with the additional thermal contributions. 

3 The Wess-Zumino model 

The Wess-Zumino model |Hj is the simplest supersymmetric quantum field theory. 
It describes a single self-interacting chiral superfield 

$ = (j) + V2 6ip + 69F 
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whose component fields are a scalar field 0, a Weyl fermion ip and an auxiliary 
scalar field F. The Lagrangian reads 

£ = *t$| wW + (W(*)L + h.c.) 

where the interaction is determined by the superpotential which we will take as 

W(Q) = -$ 2 + ^$ 3 . (8) 



Without dwelling on the superspace formalism, we rewrite the Lagrangian in 
terms of the component fields. After eliminating the auxiliary field through its 
algebraic equation of motion, the on-shell Lagrangian reads 

C = - m)tf + lid^Ad^A - m 2 A 2 ) + ^Bd^B - m 2 B 2 ) 

- gmA(A 2 + B 2 ) - \g 2 (A 2 + B 2 ) 2 - gV(A - i^B)^. 

This Lagrangian describes a scalar field A and a pseudoscalar field B, defined as 
the real components of the scalar field 0, 

A = ^(0 + 0t), 5 = --y0-0t }) (10) 

in interaction with the Majorana fermion 

* = (?> 

As a consequence of super symmetry, all fields have equal mass m, and all couplings 
are determined by the parameters g and m. We assume the coupling g to be small, 
in order to allow a perturbative treatment. 

The propagation of the three fields in a thermal background can be described 
by the usual real-time propagators for massive particles, 

Da/b(q) = — 2 it— + 2n5(q 2 - m 2 )n B (q ), 

q 2 -m 2 + ie ^ 



S(q) = (4 + m)[ j— - 27r5(q 2 - m 2 )n F (q ) , 

\q — m z + le J 

where we will denote, in a diagrammatic language, A-bosons by dashed lines and 
.B-bosons by wiggly lines, while solid lines stand for the fermion. 
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4 Super symmetry 

The action formed by the Lagrangian Q is invariant under the on-shell supersym- 
metry transformations 

5 A = IV 

SB = ^ 7 5 ^ (13) 
Sty = + m )(A + rfB)i - g(A + i-y 5 B) 2 ^ 

where £ is an infinitesimal fermionic transformation parameter. These continuous 
transformations are associated with a conserved current J M , the supercurrent. It 
is given by [TO] 

= _(0 + im)A^m + (0 - im)tfB'fV - ig{A + ij 5 B) 2 j^. (14) 
The supercharge 

« = /<!•, a*), 

lets us express the above supersymmetry transformations as the commutator (for 
bosonic operators) or anticommutator (for fermionic ones) with the charge, 

50 = -il[Q,0]±. 

In terms of the supercharge, the supersymmetry algebra reads 

{Q,Q} = 2>fP lt , (15) 

where P M is the energy-momentum operator. This relation with Poincare symmetry 
displays why supersymmetry is necessarily broken in any thermal background. The 
nonvanishing energy density of the heat bath, (Po)p ^ 0, breaks Lorentz invariance 
spontaneously and, by relation (|T5|) . also supersymmetry breaks down. 



4.1 Ward-Takahashi identities and Goldstone's theorem 

The existence of a conserved current is a highly nontrivial fact which leads to 
important relations between time-ordered correlation functions involving the sym- 
metry current. The general form of these Ward-Takahashi identities reads 

n 

&l(TJ»(x)O x { Vl ) . ..O n {y n )) = ^26^(x-y i )(T0 1 {y 1 ) . . . [Q,0i(j/i)] ± . . .O n {y n )). 

i=l 

(16) 
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It is important to note that these identities are basically operator identities, so 
that they are valid in the vacuum as well as in a thermal state, as it was es- 
tablished in f3]. Furthermore, they are valid even if supersymmetry is broken 
spontaneously. In the case of an explicit breaking, there is obviously no reason to 
expect their validity. Thus, the Ward-Takahashi identities provide a useful tool 
for the investigation of broken symmetries. 

The spontaneous breakdown of supersymmetry is characterized by some ferm- 
ionic operator O transforming inhomogeneously, 

({Q,O(y)})^0. 

The corresponding Ward-Takahashi identity, 

d?(TJ»(x)0(y)) = 5^(x - y)({Q, 0(y)}), 

can be rewritten in momentum space. By defining 

r$o(*) = Jd 4 xe^ x -y\TJ^x)0(y)), 

one obtains 

-ik^ JO (k) = ({Q,O(y)})^0. 

In order to satisfy this equation for all momenta k, the Fourier transformed cor- 
relation function Tj on the left hand side necessarily has a pole for k = 0. This 
is of course nothing but Goldstone's theorem. In the vacuum, it follows from 
Lorentz invariance that there must be a pole for all lightlike momenta k 2 = 0, that 
is, there must exist a massless Goldstone particle. One cannot however draw this 
conclusion at finite temperature. In this case, Lorentz invariance is broken and 
there can in principle be an isolated pole for vanishing momentum without the 
need for a Goldstone particle. This observation lead the authors of |U |3] to the 
conclusion that the thermal breakdown of supersymmetry is not associated with 
the existence of a Goldstone particle. A similar observation has been made for the 
related breakdown of Lorentz invariance In any case, the identification of a 
propagating Goldstone particle in a specific model requires some more effort. 

For our case of the Wess-Zumino model, the simplest Ward-Takahashi identity 
one can consider is that for a single fermion operator, 

c%(TJ»(x)*(y)) p = -i5^(x -y)m(A) p . (17) 

On tree level, this identity is clearly fulfilled in a trivial way. The behaviour in the 
interacting theory is less trivial and will be discussed in section [7| 
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Secondly, we can establish the identity for the composite mode A§> . Inserting 
the supersymmetry transformations (j!3|) into the general formula (j!6|) . we obtain 



c%(TJ»(x)A(yMz)) p = 5W(x-y)i(T*(yMz))(3 

+ 5^ 4 \x - z){$ y - im){TA{y)A{z)) p + 0{g). 

It is worth studying how this identity is satisfied at finite temperature. Let us 
switch off the interaction for a while and consider the special case y = z. In 
momentum space, the right hand side then reads 

Jd 4 xe ik ^RHS = J^L{iS(q) + (-ii-im)D A (q)) 

= -2mm J -^-^ (n F (q ) + n B (q ))5(q 2 -m 2 ). (19) 

As expected, the vacuum contributions from bosonic and fermionic propagators 
cancel because of the equality of their masses. Thus, the right hand side is trivial at 
zero temperature. At finite temperature, in contrast, the different thermal propa- 
gators for boson and fermion no longer cancel but leave a nontrivial right hand side. 
In the limits of nonrelativistic and relativistic temperatures, it is straightforward 
to calculate, 

I -2ie~ m/T ( — ^ 1 T < m 
d A xe ik{x - y) RHS = { t2 \ 2n / (20) 

—im— T ^> m. 

8 

Thereby we have neglected any corrections suppressed by factors of order T/m 
and m/T, respectively, which will always be done in the following, as long as they 
do not become relevant. 

The left hand side of (|T%|) gives 

LHS = d?(T(-$ x - im)A{x)<f*(x)A(y)*(y)) p (21) 
d 4 ^ ~-ik(x-v) f d q 



■ c 



lk{X ~ V) lAl W -H- im)D A {k - q)(-i^)S(q). 



(2tt) 4 J (2tt 

The loop integral gets two contributions from the thermal part of either of the 
propagators. The knowledge that these have support on the mass shell simplifies 
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the algebra considerably, and one calculates 



/ d'«*M MS = f** *<■¥-*->»)¥« + { . 2w)MqoW _ m >) 
J J (27r) 4 [k — qy — m l + le 

f d 4 g i{4 - m)U(U - d + m) 2 2 

f d 4 q 

= -2mm Jj^y (n ¥ (q G ) + n B (q ))5(q 2 - m 2 ). 

So, we are left with the same result as for the right hand side (fTHj) . and we have 
shown that the Ward-Takahashi identity ()18j) is fulfilled in a nontrivial way. In 
the vacuum, this would be possible only if there was a massless Goldstone particle 
generating the pole in the correlation function on the left hand side. At finite 
temperature, the situation is different. Although there is obviously no Goldstone 
fermion in the free, massive theory, nevertheless the identity is non-trivially satis- 
fied. 

The identity for the mode B^> behaves in the same way. One has 

dZ(TJ»(x)B(y)*(y)) p = d;(TJ^x)A(y)^(y)}^ 7 5 . (22) 

Comparing the relative couplings of the two modes to the supercurrent, one can 
identify the linear combination 

X = A$ - rfB^ (23) 

as the Goldstone mode. This is not surprising once rewritten in terms of the 
original chiral superfield by means of equations (jTUjl and (fTTj) . 
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{A-irfB)*=(*» ). (24) 

Thus, the Goldstone mode is the fermionic component field of the composite su- 
perfield $ 2 that determines the mass term in the superpotential (|Sj) . Note that 
the result (|20|) is nontrivial only because of this mass term. 

It will turn out in section [7| that in the full interacting theory the Ward- 
Takahashi identities are satisfied in a somewhat different way. It will be shown 
that the Goldstone mode indeed corresponds to a propagating particle. Further- 
more, the Goldstone mode x interacts with the fundamental fermion by means of 
the Yukawa interaction 

C Y = -g^(A-ij 5 B)^f. 

Thereby, a Goldstone pole will appear also in the fermion propagator which reflects 
itself in a nontrivial behaviour of the Ward-Takahashi identity (|17|). 
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So far, our discussion overlaps with earlier results from [I]. Let us go one step 
further and see what we can learn for the general case. Clearly, the way the first 
identity (|T7|) is satisfied is highly dependent on the details of the model. The second 
identity (|18|) is more universal since it involves, to leading order, only the mass and 
behaves non-trivially already at the tree level. However, we would rather expect 
a nontrivial behaviour completely independent of masses and couplings since, for 
the reasons given above, the thermal breakdown of supersymmetry is a universal 
phenomenon. According to the supersymmetry algebra (|15p. it is tied to that of 
Lorentz invariance. Consequently, the most general Ward-Takahashi identity one 
should consider is that for the supercurrent itself. Its transformation law is given 



where the additional terms vanish when taken as the expectation value in a trans- 
lation invariant state like the vacuum or a thermal equilibrium state ^Hj- This is 
basically the local form of (fTKj) as the energy-momentum tensor T^ v is the current 
generating the translations P M . The corresponding Ward-Takahashi identity then 
reads 



and is completely independent of the model. The nonvanishing of the right hand 
side is characteristic to any thermal background, and the identity implies that 
a nontrivial Goldstone mode must appear in all operators contributing to the 
supercurrent. In this way, the identity (fl"8j) is rather a special case since, in the 
massive theory, the bilinear operator A^> is part of the supercurrent (J 14)) . 

Let us verify the identity (|27)|) for our case of the Wess-Zumino model. The 
right hand side involves the thermal expectation value of the energy-momentum 
tensor. In a thermal equilibrium state, it reads 



by 



{Q,J»} = 2T^ lu + ..., 



(25) 



d;{TJ^x)r(y)) = S^(x - y) 2(T^) %1 



(26) 



(Tn/3 = dmg(p,p,p,p), 



where the energy density p is given by 




) 



3/2 



d 3 g 



4me -m/T 



P 



(2tt) 3 



E q (2n F {E q ) + 2n B (E q )) 




T»m. 



The pressure p is related with the energy density as 
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Hence, on the right hand side of the identity (|26|) we find, for T ^ m, 



2(2**% 7 „ = 8me-/ T (^X^ { ^ ) . (27) 



2^ 



T 

— 7 



For T>m, we have 

2(T^),7, = T T 4 ( /I ]. (28) 

This must be compared with the left hand side of (|2fij). 

LiLS = d£(T(-$ x - im)A(a;)7^(a;)tf Q/M-ft, + im)A(y)) p 

+ ^(T(^ - im)i 1 h B{x) 1 ^(x)^(y)Yil" , ^ y + im)B{y)) p . 

It can be evaluated in the same manner as equation ([21)1 . The same cancellations 
take place, and one obtains, after some algebra, 



J d*xe lk ^LHS = 4 j-0- q ^27r(n F (q o ) + n B (q ))6(q 2 - m 2 ). 

Calculating the remaining integral, one ends up with exactly the same result as 
for the right hand side, equations (|2*7j) and (|28jl. Thus, the Ward-Takahashi iden- 
tity (|2fjjl is shown to be fulfilled in a nontrivial way. 



4.2 Supersymmetric sound 

In view of the fact that the quantum field theoretical framework not necessar- 
ily requires a propagating Goldstone particle associated with the thermal break- 
down of supersymmetry, it is even more interesting that the existence of a mass- 
less fermionic excitation can be deduced in a complementary setting. In [S| (see 
also |12j). a hydrodynamic approach to the supersymmetric plasma has lead to 
the prediction of a slow-moving collective excitation whose existence should be as 
general as that of sound waves. 

Consider a relativistic perfect fluid with an energy-momentum tensor of the 
form 

(TH/3 = diag(p,p,p,p). 

If the system is disturbed by a small, spacetime-dependent variation of, say, 
the temperature, AT(x), the conservation of the energy-momentum tensor, 
d fl (5T' J ' u )i3 = 0, translates to a wave equation for the disturbance, 

(|^ 2 - |^ 2 ) AT(x) = 0. (29) 
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It describes the propagation of sound waves with the velocity 



Thus, as a consequence of the breakdown of Lorentz symmetry in the thermal bath, 
small perturbations in the energy-momentum tensor propagate as sound waves. 
Their quanta, the phonons, can be viewed as the Goldstone bosons associated 
with the spontaneous symmetry breaking. 

Now, in the supersymmetric case, one can imagine the system undergoing a 
small, spacetime-dependent supersymmetry variation £(x). The conservation of 
the supercurrent, d^{SJ^(x))/3 = 0, translates to a wave equation for the transfor- 
mation parameter, 

(p7odo+Pjd)Z(x) = 0, (31) 

where the components of the energy-momentum tensor enter through the trans- 
formation law (|25j) . This Dirac equation describes a massless fermionic excitation 
propagating with the velocity 

vss = -• (32) 

Thus, in a system whose supersymmetry is broken by the thermal bath, there 
should exist 'supersymmetric sound waves'. Their quanta, interpreted as the Gold- 
stone fermions associated with the spontaneous symmetry breaking, are naturally 
called phoninos. 

Both sound and supersymmetric sound have a very characteristic dispersion 
law given by 

f - T«m 

v 2 s = vss={ r T „ ( 33 ) 

in the non-relativistic and relativistic limits. 

Despite the similarity to sound waves, there is surely no classical picture of 
these fermionic waves. It should however be possible to verify their existence 
in the framework of thermal field theory. Just as sound waves appear as poles 
in the correlator (T^^T^y)) p, supersymmetric sound must lead to poles in 
the correlation function (J tt (x)J„(y))p. Since supersymmetric sound appears as a 
collective phenomenon, it is expected to appear in the quantum field theoretical 
framework only as a nonperturbative effect. 



5 One-loop corrections 

Though the thermal breakdown of supersymmetry is evident already at the tree 
level, a deeper understanding certainly requires the investigation of higher-order 



12 



effects. The one-loop effects have already been investigated in earlier works, so 
that we can refer to f3J |Hj for explicit calculations. 

Firstly, the interaction with the heat bath leads to the appearance of a non- 
trivial expectation value of the fundamental scalar field. Its value is given by 

{—gma T <^ m, 
gT 2 _ ^ (34) 
"2^ T>>m ' 
where we have introduced the small parameter 

-^/^^-^sr <35> 

The nonzero value is another clear sign for the breakdown of supersymmetry. It 
must be stressed, however, that it is not a necessary condition like in the vacuum. 
In the massless Wess-Zumino model, for example, the expectation value vanishes 
because of chiral symmetry while supersymmetry is nevertheless broken for the 
general reasons given above. 

Secondly, the interaction leads to a modification of the propagators. The dom- 
inant self energy corrections come from the proper one-loop diagrams and, with 
a contribution of the same order of magnitude, from the coupling to the thermal 
expectation value of the scalar field. For momenta p = (m, 0), the relevant real 
parts of the thermal self energy functions Ilf (p) have been calculated in Their 
result can easily be generalized to momenta on the mass shell, since none of the 
diagrams shows a dependence on the three-momentum p, as long as p 2 = m 2 . As 
a consequence, the self energy corrections only lead to a small displacement of the 
masses, 

m 2 = m 2 + nf(p)| p2=m2 
without changing the dispersion relation 

p 2 Q = p 2 + m 2 . 

Explicitly, the values of the effective thermal masses are given by 

r m \ = m 2 -2g 2 T 2 

T <^ m : { m% = m 2 T ^> m : < m\ = m 2 

I m% = m 2 — g 2 T 2 

(36) 

As yet another evidence for the breakdown of supersymmetry, we find that the 
mass degeneracy of the three fields is lifted by the one-loop thermal self energy 



m\ 


2 

= m 


- ^g 2 am 2 


m% 


2 

= m 




m% 


2 

= m 


— 2g 2 am 2 
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Figure 1: One-loop contributions to the fermion self energy 



corrections. As expected, the mass splitting vanishes in the limit of zero temper- 
ature when supersymmetry is restored. 

In a model where supersymmetry is broken spontaneously already in the vac- 
uum, the masses of the members of one supermultiplet are related by the mass 
formula 

^(-1)^(27 + l)m»=0. 
j 

For the effective thermal masses (|36|). we find that this relation is no longer valid, 
at least for low temperatures where m 2 A + m 2 B 7^ 2m\. This raises the question 
whether the breakdown of supersymmetry should at all be considered spontaneous 
or rather explicit. 

6 The phonino 

In the previous section, we have seen that the influence of the heat bath on the 
fundamental fields is basically a small shift in their masses. According to the 
general arguments given in section 0J however, we suspect the existence of addi- 
tional fermionic excitations. Hence, we should take a closer look at the full self 
energy function that determines the poles of the full fermion propagator. These 
are characterized by the condition 

det S'\p) = -idetty-m-Eptp)) =0. (37) 

To leading order, the self energy gets contributions from the coupling to the 
expectation value of the scalar field A as well as from the proper one-loop diagrams 
as drawn in figure d Altogether, one finds 

-i E p (p) = -2ig(A) fi - V J-j0- D(p - q) (S(q) + vfS(q)vf) . (38) 

We are particularly interested in the limit of small momenta. Thus, we set the 
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three-momentum to zero and calculate 



E(„ Q)-2a(A) + 8jV f A3q ( n *(E,) 2E* t -ri ny (E q ) 2£? 
S M -2 g (A), + — J— y ^ , + - , 



where E q = \/q 2 + m 2 . One observes that the second part develops a pole for 
vanishing momentum. The reason is of course the degeneracy of bosonic and 
fermionic masses. In the limit p — > 0, the poles of both propagators in the loop 
coincide, leading to an anomalously big self energy. 

In the relevant limits of low and high temperatures, the integration can easily 
be performed. For small momenta po, the result can be expressed in terms of (A)p 
as 

E,(p ,0) = 2 9 (^-^V. 

Po 

Inserting this result into condition (|H7|) for the propagator poles, one obtains the 
dispersive equation 



with the solutions 

\p \=m + g(A)p + O(g 2 (A} 2 ) 1 
\p \=g(A), + O(g 2 (A)l). 

The first solution reproduces the small shift in the fermion mass that was already 
calculated in the previous section. The existence of a second solution indicates 
another excitation with a tiny mass. In fact, such a pole is what we expect, albeit 
with exactly vanishing mass. It thus seems that our simple one-loop calculation is 
not sufficient. 

As it will turn out in the following sections, a consistent calculation of the 
full propagator requires to take into account the full propagators in the internal 
lines, that is, to consider the differences in the effective masses of bosons and 
fermion due to the breakdown of supersymmetry. In this way, one must perform a 
nonperturbative, selfconsistent calculation in order to uncover the desired massless 
pole in the fermion propagator. 

6.1 The case of high temperature 

Let us now analyze the fermion self energy more carefully. At first, we restrict 
ourselves to the relativistic case. By the assumption that the mass shifts are 
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small compared to the mass, our discussion is then limited to a temperature range 
m <T <C g~ x m. 

Since the interaction with the heat bath modifies the propagators by a small, 
constant mass correction, we will approximate the full propagators by the one-loop 
resummed propagators 



S(q) = {4 + m + E'(q)) (- \ r — - 2n5(q 2 - m|)n F (g ) J 



(39) 



^A/B(q) = — 2 2 r~ + 27T6 (l ~ m A/B) n B(qo 

q z - m\ /B + ie ' 



with the effective thermal masses as given in equation (|36|) . As one can check 
from equation (|5Bj). the one-loop self energy splits into a constant part 
proportional to the unit matrix and a traceless part E 7 (which is an odd function 
of the momentum), so that one has 

S'(g) = S 1 -S 7 (g). 

With this approximation for the full propagators, we can now evaluate the 
fermion self energy ()38jl . For small momenta k, the product of both propagators 
in the loop can be approximated as 



A(Jfc - q)S(q) = -2m (4 + m + E'fa)) 5(g 2 - m 2 ) " B 2 w; J 2 " F ^ , (40) 



| ng(go) + ^f(<?o) 
m\ — mf — 2qk ' 



where we only keep corrections linear in k to small terms of order g 2 . A general 
loop integral then reads, for sufficiently regular /, 

<V<I V l {k-q)S{q)f{q) 



(2nY 

_. f d 3 g ^ n B (E q )+n F (E q ) j + m + Xjq) 

—+f 2E i ml-m 2 -2qk nq) - [ * L) 

One finds that the small mass splitting between bosons and fermion regularizes the 
loop integral. So, in contrast to the one-loop calculation with tree-level masses, 
the integral no longer diverges in the limit k = but becomes anomalously large 
because of the small mass difference appearing in the denominator. We assume 
that qk m\ — mf with a typical loop momentum q ~ T which restricts ourselves 
to momenta k -C g 2 T. 
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The explicit calculation of the full self energy (J55J) then gives 



»4-mV '7(2>r) 3 E 



4# 2 / d 3 g n B (£ 9 ) +n F (£ g ) 



8# 2 (&q7o + 



d 3 g 



(2*) ! 



E q (n B (E q )+n F (E q )) 



One can now insert the mass splittings resulting from equation . 

2 2 2 2 2ti2 

my - m A = m B - = g T , 
and perform the integration. One ends up with the simple result 



Ep(k) = -m - ^— (ko'jo + §£7) 



Remarkably, the result is exactly the negative Lagrangian mass so that, in our 
approximation, the full inverse fermion propagator reads 



2 

iS-\k) = ff-m- Zf,(k) « ^ (fc o7 o + |fcf) • (42) 



Because of the cancellation of the mass term, we find that the full fermion propa- 
gator has an additional pole for vanishing momentum as well as on the dispersion 
curve 

ko = ±l\k\. (43) 

This displays the existence of a massless excitation propagating with velocity 
v = |, just as one expects from the general arguments given in section l4~2l Con- 
sequently, we identify this pole in the fermion propagator as the phonino. 

It should be added that this mechanism to cancel the fermion mass by loop 
corrections goes back to a widely unnoticed paper by Kapusta [H] who performed 
a similar calculation for the somewhat simpler case of the massless Wess-Zumino 
model. However, the author attributed the existence of the soft fermion to chiral 
symmetry and not to the breakdown of supersymmetry. 
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a 



Figure 2: Tree-level Yukawa coupling and one-loop corrections 
6.2 The case of low temperature 

The above calculation can in principle be translated also to the case of low temper- 
ature. It turns out, however, that this case requires a more thorough investigation. 
We have seen that contributions of the form (|40j) from two propagators of nearly 
degenerate mass and opposite momentum lead to an anomalously large value of 
the self energy. So, for consistency, one must take care of all diagrams of a similar 
structure that could give equally large contributions. 

Take a look at the vertex corrections in figure 121 For vanishing external mo- 
mentum and low temperature, the contribution from diagram b is of order 

3 f d 3 q n F (E q ) + n B (E q ) m I 



(2tt) 3 E q m% — m 2 A m? ® 

Hence, it is as relevant as the tree-level coupling. For high temperature, in contrast, 
one finds 

3 f d 3 q n F {E q ) + n B {E q ) T 1 m 
a m I - — — - - * n — n ~ — Q, 

so it does not give a relevant contribution because of the tree-level Yukawa cou- 
pling being proportional to m C T. The same is true for diagrams with fermion 
exchange like in figure E]c. Thus, for high temperature there is no need to take 
these diagrams into account, while for low temperatures one must evaluate all 
relevant contributions to the vertex. 

Obviously, many-loop diagrams play an equally important role, as long as each 
intermediate vertex or line is again followed by a pair of particles with a tiny mass 
difference, which leads to a cancellation of the small coupling constants and dis- 
tribution functions that otherwise suppress any higher-order diagram. Altogether, 
one must sum up all ladder diagrams of this kind. 

The standard way to perform such a resummation is by solving the correspond- 
ing Bethe-Salpeter equations. Their general form reads, diagrammatically, 
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Here, the shaded spot stands for the full proper (one-particle irreducible) vertex 
while the square denotes any proper (two-particle irreducible) scattering subdia- 
gram, and a sum over all possible intermediate states is implied. The internal lines 
stand for full propagators. 

In our specific case, the intermediate state in the Bethe-Salpeter equation is al- 
ways a boson-fermion pair of nearly degenerate mass. The two-particle irreducible 
scattering amplitude gets contributions from both boson and fermion exchange. 
So, the Bethe-Salpeter equation for the vertex function involving the A boson reads 





(44) 



A similar equation must hold for the B boson vertex, 





(45) 



Let us translate these diagrammatic equations into mathematical language. 
The amputated full proper vertex functions will be denoted by T^kjp) and T B (k,p), 
respectively, 




= T A (k,p), 




= r B (k,p), 



which we want to calculate for small momenta k. The thermal contributions to the 
Bethe-Salpeter equations come from loop momenta on the mass shell. Therefore, 
we consider only on-shell momenta p 2 = m 2 . 

All contributions to the Bethe-Salpeter equations have a very similar structure. 
Approximating the full propagators in the internal lines by the one-loop resummed 
propagators we again find expressions of the form (|4*Tj). For simplicity, we 

introduce 



A ^4 = niz, — m 



i^ 
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2 2 

g am , 



A 



13 



2 2 o 2 2 



(46) 
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Like in the case of high temperatures, we only consider momenta k small compared 
to the mass differences and keep only the linear terms, while small corrections 
of order Aj/m or k/m will be neglected. Furthermore, we can safely neglect 
contributions of order pjm or q/m (where q is the momentum of the fermion 
in the loop), since p and q are thermal momenta of order y/Tm which, for low 
temperature, is negligible compared to the mass. The propagators of the exchanged 
particles are then essentially constant but depend on the relative sign of the zero 
components of outgoing and loop momenta. 

Altogether, the Bethe-Salpeter equation for the A boson vertex can be written 

as 



r A (k,p) = -2ig-4g 2 



d 3 q e~ Eq / T 



(2tt)3 E q 



+ 



if) + 2m 3m 
3m 2 — m 2 

+ 2m 3m 

T- + 



-wr 



$ + 171 

A A - 2qk 

q 1 + m 



3m 2 / A A — 2qk 



F A (k,q) 
F A (k,q) 



q =E q -signp 



q =-E q -signp 



—f 1 + 2m m 
3m 2 — m- 

—p 1 + 2m m 



—q + m 



A 



B 



2qk 



ij 5 r B (k,q) 



-fi + m . 5 



— m Zj 3m 2 / Ab — 2qk 
For the B boson vertex, one finds 

d 3 q e~ Eq/T 



irr B (k, q ) 



q =E q -signp 



q =-E q -sign p 



i-fr B (k,p) = -2ig-4g 2 



3rri 



+ 



(2tt) 3 E, 
-m 



+ 



q 

$ + m 
-m- I A \ — 2qk 

in \ $ + m 



3m 2 



+ 



3m 2 / A^ — 2qk 
m 



F A (k,q) 
T A (k,q) 



qo=E q -signp 



qo=-E q -sigapo 



+ 



in 



+ " l r^Ak.q) 



-171" 



3m 2 7 A 



A B - 2qk 
—d + m 



2qk 



h 5 r B (k, q ) 



qo=E q -sign p 



q =-E q -signp 



(47) 



(48) 



These seemingly complex coupled integral equations can be simplified consid- 
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erably by introducing the functions 

Rewriting the Bethe-Salpeter equations in terms of these functions, it turns out 
that the right hand sides involve only the symmetric parts 

VaA^p) = (t A , B (k,p)+f AiB (k,-p)) . (49) 

V / po=E p 

By straightforward calculation, one can express the Bethe-Salpeter equations in 
terms of the functions Va,b as 

(50) 

where we have again made a linear approximation for small momenta k. 

For vanishing k, one finds that the functions Va, B {0,p) do not depend on the 
momentum p at all. The integration can then easily be performed by making use 
of the relations 

V f d 3 q _ Eq/T 1 2 3A, 1 A fl 

/ 6 = — G OLTYl = = 

m 2 y (2tt) 3 2 y 28 m 4 m 

which immediately follow from equations (|35J) and f)46|) . One arrives at the system 
of linear equations 



(51) 



Va(0, P ) = ^ \-2ig + ^— (2V A (0,p) + \v B %p) 
Aa { 28 m \ 3 

V B (0,p) = -— -2ig - --^ --Va (0,p + -V B (0,p) 
A B { 4 m \ 3 3 

which is easy to solve. The solution is given by 

V A (p,p) = V B (p,p) = 1 r--2ig. (52) 

It is illustrative to express this result in terms of the amputated vertex functions 
Ta and T B we started with. One finds 



r B (0,p) = -2ig (-i 7 5 ) , r A (0,p) = -2ig ("^£) • 



(53) 
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In other words, the Yukawa coupling to the B boson comes out of the whole 
resummation without a change while the coupling to the A boson is enhanced by 
a factor of 

_Aa = 7 
A B 3' 

Let us now calculate the momentum dependence of the vertex functions. To 
this end, we expand in k, 

V A , B (k,p) = V AB (0, P ) + k^M, 

which is substituted in the Bethe-Salpeter equations (J5UJ) . The resulting integral 
equations for the functions C^ B (p) can be simplified by replacing the known values 
of Va ; b(0,p) and integrating, so that one obtains 

" " & /<& ^ ("I ^ + 1 H ■ 

(54) 

At this point, it is no longer possible to neglect the momentum dependence since in 
particular the vector components Ca,b(p) receive an important contribution rising 
quadratically with \p\. Yet, the system can be solved by multiplication by e~ Ep ^ T , 
followed by integration. In this way, one finds a linear system with the solution 

/|K' /T VM - f^- E -' T ^M - ^ (m*» + T*7) . 

(55) 

One can now proceed and calculate the functions C^ B (q) by inserting this result 
into equations (|53j). We skip this last step since the result will not be needed in 
the following. 

Instead, we proceed to the calculation of the full fermion propagator. With our 
knowledge of the full one-particle irreducible vertex functions, we can now evaluate 
the full proper self energy. In terms of full propagators and vertex functions, the 
main contribution from ()38|) translates to 

-i V fi (k) = ~2ig J-0- 4 V(k- q) (T A (k, q)S(q) + T B (k, q)S(q)z 1 5 ) , 

which, by making use of equations (}4*Tj) and (|49jl. is easily brought into the form 

Z fi (k) = -2ig j^- 3 e -^(V A (k,q) + V B (k,q)). 
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With our results for Va,b in (jo^j) and (jo3j) . one calculates 

Ep(k) = -m 1~ ( A; 7o + — ^7 J • (56) 

gra \ my 

Thus, the self energy £^(0) = — w exactly cancels the Lagrangian mass term in 
the full inverse propagator which thereby becomes 

iS-\k) =#-m- E p (k) « ( fc o7 o + • (57) 

So we have finally found the additional massless pole in the full fermion propagator. 
It becomes singular for vanishing momentum as well as on the dispersion curve 

k = ±— \k\ (58) 
m 

which is nothing but the dispersion law predicted in section 14.21 Therefore, this 
pole can be identified with the phonino. Thus we have shown that the phonino pole 
in the fermion propagator is present at any moderate temperature, only the residue 
and dispersion law change with temperature. In the limit T — > 0, the residue 
vanishes and the Goldstone mode disappears, since supersymmetry is restored. 

Due to our approximation, we have proven the existence of the phonino pole 
only for small momenta. We have assumed that kq <C with a typical ther- 
mal momentum \q\ ~ yTm, so that the derivation is restricted to momenta 
|fc| <C A a/ 'yTm. Interestingly, this is not only a technical point but of physical 
relevance. In the hydrodynamic picture, the existence of sound and supersymmet- 
ric sound waves is restricted to long wavelengths much greater than the mean free 
path so that there is always time to establish local thermodynamic equilibrium. 
At higher frequencies, the waves are strongly damped as it was shown in 

The fact that our calculation requires a resummation of both one-loop correc- 
tions to the propagators and higher-order vertex corrections shows that the exis- 
tence of the phonino is really a nonperturbative phenomenon. This is consistent 
with the interpretation as a collective excitation according to the hydrodynamic 
explanation given in section l4~2l Furthermore, it explains why earlier calculations 
were not able to relate the Goldstone mode to a propagating excitation by simpler 
one-loop calculations [HE]- 

The existence of the low-temperature phonino pole and dispersion law was 
already shown in [H] by investigating the homogeneous Bethe-Salpeter equations for 
the full (i.e., not one-particle irreducible) vertex functions. Our explicit calculation 
of the full fermion propagator is however a new result which will be essential for 
the following discussion of the Ward-Takahashi identities. 
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7 Verification of the Ward-Takahashi identities 



The appearance of the massless phonino derived in the previous section is appar- 
ently linked with the breakdown of supersymmetry. In order to obtain a complete 
picture and to prove that the phonino is indeed the Goldstone fermion of sponta- 
neously broken supersymmetry, let us now take a closer look at the supersymmetric 
Ward-Takahashi identities and verify that the phonino contributes in exactly the 
way expected for a Goldstone particle. 

In section 01 we have derived the identity 

d*(TJ»(x)% (y)) p = -im(A) p 8®{x - y). (59) 

At the one-loop level, the right hand side is nontrivial because of the nonvanishing 
thermal expectation value of the scalar field given in equation . Thus, the 
left hand side must get a contribution from a Goldstone mode. Surely, we expect 
the phonino to act this part, with the main contributions coming from the class of 
diagrams drawn in figure H3 Here, the encircled cross denotes the supercurrent (1141) . 
The shaded spot stands for the full vertex, and the thick lines mean full resummed 
thermal propagators. 

The evaluation of the these diagrams thus gives, to leading order, 

(TJ^x)^(y)) p = 

I J0f e ~ lk(x ~ y) {/(& ~ H ~ im) YVA{k ~ ^ s ^ Ta ^ 

4 I (60) 

Going over to momentum space, we amputate the fermion propagator and 
define 

T%(k) = Jd*xe*^(TJ''(x)%(y))f l S- 1 (k) 




\ / 

N X 




Figure 3: Main contributions to the Ward-Takahashi identity (|59j) 
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so that the Ward-Takahashi identity (|59|) we want to verify can be written as 

-i* M Ijff(*0 = -im{A)pS-\k). (61) 

Now to the calculation of the left hand side. The loop integrals in (|5Uj) are 
again of the form (j41j) well-known by now, and the Dirac algebra is simplified by 
the fact that the loop-momentum is on-shell. Since we are mainly interested in 
the role of the Goldstone mode, we consider small momenta k and only keep the 
linear terms. We find 

d3 9 \ i oj. 7 sn B (E q )+n F (E q ) T A {0,q) 



-ik,T%(k) = -i f±jL (-2iqk) 



2E g m% - m\ 



+ 24 qk 



n B (E q ) + n F (E q ) tj 5 T B {0,q) 
2E q ml - m% 



The integration gives, for nonrelativistic temperatures T ^ m, 

-ik,T%{k) = — k olo + -kj). (62) 
9 V m J 

Similarly, for the relativistic case T»m, one finds 



7T 



2 T 2 / j 



-ik l ^%ik) = ^fco7o + |^7j • (63) 

One immediately observes the similarity to the inverse phonino propagators ob- 
tained in equations (jl2| and (jSZJ), and indeed, one finds by comparison with equa- 
tion ()34j) that the difference is just a factor of —im{A)p. So, the Ward-Takahashi 
identity (|6T|) is fulfilled through the contribution from the phonino, and we can 
conclude that the phonino is the Goldstone fermion associated with the sponta- 
neous breakdown of supersymmetry in the thermal background. 

Once we have convinced ourselves that the Ward-Takahashi identity is satisfied 
at finite temperature in the full resummed interacting theory, we can in turn use it 
to obtain more information about the full propagator. If one keeps also the linear 
terms in the calculation of Tj^(k), equation (|BDj). one can deduce the nonlinear 
corrections to the phonino dispersion law. In agreement with [6 , they are found 
to be small in the region where the phonino exists. 

Let us now proceed to the second Ward-Takahashi identity, the one for the 
composite mode that was derived in equation (fT%|) . Its right hand side was 
evaluated already in equation (|20|) so that we can rewrite the identity in terms of 
(A)p as 

aZ(TJ»{z)A{y)V {y)) fi = i^(A) fi S^(x - y). (64) 
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Again, the right hand side is nontrivial at nonvanishing temperature, so that the 
left hand side must get a contribution from the Goldstone mode. 

The diagrams contributing to the left hand side fall into two classes, drawn 
in figure EJ For the first class of diagrams, one must have control over the full 
one-particle irreducible four-point function denoted by a shaded box that could 
in principle have additional poles. However, in the limit of small momentum we 
are interested in, it can easily be shown that this is not the case. The relevant 
contributions to this four-point function can be resummed by means of the same 
Bethe-Salpeter equations that were set up for the calculation of the one-particle 
irreducible three-point functions in section lfi~2l At this point, it is only important 
that this resummation does not introduce any new poles, in other words that the 
homogeneous Bethe-Salpeter equation does not have any nontrivial solutions. This 
can easily be seen from equation (|51|) which, as a homogeneous system, has only 
the trivial solution. Hence, this class of diagrams stays finite in the limit of low 
momentum. 

Thus, the relevant contributions come from the second class of diagrams in- 
volving the phonino, schematically drawn in figure Ejb. Here, the big smudge 
stands for the sum of all diagrams discussed above in connection with the first 
Ward-Takahashi identity. So, for the evaluation of this correlation function in the 
low-momentum limit, we can take advantage of our earlier results and write 



where the numerical factor is the same for low and high temperatures. 

We can now make use of the first Ward-Takahashi identity (|6*T|) and obtain 





JA*(k)= i —(A) f} 



(65) 




a 



b 



Figure 4: Contributions to the Ward-Takahashi identity (j64j) 
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which is nothing but the Ward-Takahashi identity (JMj) we intended to prove, 
written in momentum space. 

Let us compare the ways the Ward-Takahashi identity is fulfilled for the free 
and interacting theories. In the free theory, the loop diagram in figure E]a (without 
any interaction) is the only one that contributes to the left hand side of the identity. 
Because of the equality of the masses of both propagators, it becomes singular in 
the vanishing momentum limit, thereby saturating the Ward-Takahashi identity 
non-trivially without the need for a massless mode. In the interacting theory, 
in contrast, there appears a mass splitting between fermion and bosons, so that 
diagram 0]a stays finite even in the limit of vanishing momentum. Hence, there 
must be a massless Goldstone mode that contributes to the left hand side in order 
to saturate the Ward-Takahashi identity. Our calculation shows that the phonino 
pole in the fermion propagator exactly does this job. 

Finally, the Ward-Takahashi identity involving the supercurrent, 

d?(TJ»(x)J»(y))p = 6®(x - y) 2<T"% 7/i , (66) 

can be verified in the same way. The dominant contributions to the left hand side 
come again from diagrams with a single fermion intermediate state. Hence, we can 
write, in momentum space, 

rfj(k) = Jd 4 xe ik ^(TJ»(x)r(y)) p 

= r%(k)S(k)T^(k). 

By making use of the Ward-Takahashi identity (}6*T]) . we find 

-W) = -im(A)pf^(k), 

which can now be calculated by using our earlier results on Tj^(fc). 

First, in the limit of low temperature, we make use of (Jo^j) and obtain the left 
hand side of the Ward-Takahashi identity as 

- iVTj(k) = -^>|g ( £ ) = Sme-^ ( 7o_ ) 

For high temperatures, one calculates in the same way, using (jHHJ), 

Comparing with our earlier results on the energy-momentum tensor, equations (}2Tj) 
and (jZEj) , we find that this can be written as 

-ik,Yfj{k) = 2{T^) plix . (67) 

This is however nothing but the momentum space version of the identity (JBTIj) we 
wanted to prove. 
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8 Conclusion 



Our calculations have lead to a consistent picture of the behaviour of the Wess- 
Zumino model at finite temperature. At moderate temperatures, the interaction 
with the thermal background leads, on the one-loop level, to a small splitting of 
the effective masses, and the scalar field develops a nontrivial thermal expectation 
value. Both are signs, if not necessary conditions, for the breakdown of supersym- 
metry. As an additional feature, the phonino pole appears in the propagator of 
the fundamental fermion, as a full nonperturbative calculation shows. The rea- 
son is the continuous interaction with boson-fermion pairs of nearly degenerate 
mass. As a consequence, the same pole appears in all modes that couple to the 
supercurrent. This relation to the supercurrent indicates the role of the phonino 
as the Goldstone particle of spontaneously broken supersymmetry, which could be 
confirmed by the investigation of the Ward-Takahashi identities that are saturated 
by the contribution from the phonino in a nontrivial way. 

However, our results are not restricted to the Wess-Zumino model. The general 
picture we have drawn and verified in this simple model, according to which the 
breakdown of supersymmetry is a consequence of the breakdown of Lorentz invari- 
ance, can immediately be translated to any supersymmetric model. Just as the 
existence of supersymmetric sound was derived in jB] in a model-independent way, 
the Ward-Takahashi identity for the supercurrent generically predicts the existence 
of a Goldstone pole in any mode that couples to the supercurrent, as long as the 
fields involved contribute to the nonvanishing energy density of the thermal bath. 
By the same mechanism that we explicitly studied for the Wess-Zumino model, 
this pole must be associated with a propagating particle whenever the interaction 
lifts the boson-fermion mass degeneracy that otherwise allows the saturation of 
the Ward-Takahashi identities without the need for a Goldstone particle. 

Up to now, our analysis was restricted to chiral superfields and (in the vacuum) 
unbroken supersymmetry. In order to come closer to phenomenologically attractive 
models, it is necessary to understand the behaviour of gauge fields [TH] as well as 
the interplay with other mechanisms of supersymmetry breaking. 
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